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The Josephson effect is generally described as Cooper pair tunneling, but it can also be understood 
in a more general context. The DC Josephson effect is the pseudo-Goldstone boson of two coupled 
systems with a broken continuous abelian (7(1) symmetry. Hence, an analog should exist for systems 
with broken continuous non- Abelian symmetries. To exhibit the generality of the phenomenon and 
make predictions from a realistic model, we study tunneling between antiferromagnets and also 
between antiferromagnetic d-wave superconductors. Performing a calculation analogous to that of 
Ambegaokar and Baratoff for the Josephson junction, we find an equilibrium current of the staggered 
magnetization through the junction that, in antiferromagnets, is proportional to §_l x s_r where sl 
and s_r are the Neel vectors on either sides of the junction. Microscopically, this effect exists because 
of the coherent tunneling of spin-one particle-hole pairs. In the presence of a magnetic field which is 
different on either sides of the junction, we find an analog of the AC Josephson effect where the angle 
between Neel vectors depends on time. In the case of antiferromagnetic d-wave superconductors we 
predict that there is a contribution to the critical current that depends on the antiferromagnetic 
order and a contribution to the spin-critical current that depends on superconducting order. The 
latter contributions come from tunneling of the triplet Cooper pair that is necessarily present in 
the ground state of an antiferromagnetic d-wave superconductor. All these effects appear to leading 
order in the square of the tunneling matrix elements. 



I. INTRODUCTION 

One of the most striking manifestations of supercon- 
ductivity is the Josephson effect, which consists in the co- 
herent tunneling of Cooper pairs across a junction when 
both sides are superconductors. Interestingly, this oc- 
curs even though the Hamiltonian contains only single- 
electron tunneling. The resulting current is proportional 
to sin((^) , where (p is the difference between the phases of 
the two superconducting order parameters. In the pres- 
ence of a constant potential difference, gauge invariance 
requires the phase difference to increase linearly with 
time, leading to an alternating current whose frequency 
depends only on universal physical constants and not on 
any material parameter. 

Superconductivity is just one example of coherence and 
spontaneous symmetry breaking. The order parameter is 
a complex number so the broken symmetry is U(l). The 
superconducting state selects a phase and the Josephson 
effect arises from the tendency to make the phase uniform 
across the tunnel junction. 

Since there are many other types of order and corre- 
sponding broken symmetries, the question of the analog 
of the Josephson effect in such cases arises naturally^.. 
Indeed, one should expect differences in order parame- 
ters across a junction to lead to a coherent tunneling of 
the condensed objects that exist in the broken symmetry 
state. This question is especially relevant in the context 
where junctions between magnetic materials are of ut- 
most importance for spintronics. As a matter of fact, 
theoretical predictions have been made concerning the 
possible existence of equilibrium spin currents^ of the 
Josephson type in ferromagnetic (FM) tunnel junctions, 



in analogy with superconducting junctions— ~—. FM long 
range order is a realization of spontaneous 50(3) sym- 
metry breaking and an equilibrium spin current would 
result from the exchange coupling between the magnetic 
moments in the two leads, which favors alignment of or- 
der parameters. From a Ginzburg-Landau point of view, 
there is a term proportional to Mi - M^j where Ml and 
Mr are, respectively, the magnetic moments on the left 
and on the right of the junction. The Heisenberg equa- 
tions of motion thus lead to ~ Ml X Mr, corre- 
sponding to a spin current. 

In this paper, we study Josephson-like phenomena be- 
tween antiferromagnets (AF), normal and superconduct- 
ing. Such phenomena are interesting for several reasons. 
In the context of spintronics, it has been shown experi- 
mentally that the absence of a net angular momentum in 
AF results in orders of magnitude faster spin dynamics 
than in FM, which could expand the now limited set of 
applications for AF materials^. More generally, in AF the 
Neel order parameter breaks both lattice translation sym- 
metry and SO(3) spin rotation symmetry so the situation 
is less straightforward than in FM. In addition, AF are of- 
ten close to superconducting phases, as is found in heavy 
fermions, high-temperature superconductors and layered 
organic superconductors. Understanding Josephson-likc 
phenomena in AF is the first step towards more general 
studies with coexisting antiferromagnetic and supercon- 
ducting order parameters, which we also pursue in this 
paper. Generalized Josephson effects may help identify 
homogeneous coexistence of antiferromagnetism and su- 
perconductivity in real materials. Studies along these 
lines have recently appeared for ferromagnetism coexist- 
ing with p— wave superconductivity^. 
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We begin in Sec. II from a microscopic mean-field 
model for itinerant antiferromagnets and for single elec- 
tron tunneling and establish analogies with the BCS 
ground state. This microscopic calculation leads in Sec. 
Ill to an explicit expression for the analog of the critical 
current and its temperature dependence. In addition, we 
show that (a) Cooper-pair tunneling is replaced by tun- 
neling of a spin-one neutral particle-hole pair (b) time de- 
pendence introduced by external magnetic fields resemble 
the AC Josephson effect but there are many differences 
because of the non-Abelian nature of the problem and be- 
cause spins do not couple to the gauge field but directly 
to the magnetic field. We also discuss in Sec. IV how to 
observe this effect by coupling to a ferromagnet through 
a tunnel junction. We then move in Sec. V to the case 
of d— wave superconducting antiferromagnets and show 
that there is a contribution to the Josephson charge cur- 
rent that is modulated by the antiferromagnetic order 
parameter and, conversely, a contribution to the spin- 
Josephson current that is modulated by the supercon- 
ducting order parameter. Our results are summarized in 
the discussion and conclusion section. Additional details 
of the calculation may be found in Ref^. For higher- 
order effects in the tunneling matrix element s 10 ! 11 that 
we do not discuss here, see also Ref. H 



II. MODEL 

A. Antiferromagnetic state and analogy with BCS 
ground state 

The Hamiltonian for a tunneling junction consisting of 
two leads of an AF material and an insulating barrier 
between them reads 



(1) 



where -ffz,(i?) is the Hamiltonian of the left (right) AF, 
Ht is the tunneling part connecting the two leads, and 
c ko-( Ckcr ) an d ^qo-^qc) are the fermion creation (annihi- 
lation) operators of the left and right leads, respectively. 
In the following discussion the quantum numbers k and 
q will also denote implicitly the left and right lead. 

We model the AF on each side of the junction by a 
one-band Hubbard Hamiltonian treated in the Hartree- 
Fock approximation for a static spin-density wave (SDW) 
with wave- vector Q. Without loss of generality, we as- 
sume the SDW mean field to be polarized along the spin 
quantization axis. Following Ref. 12 --, we write 



Hr 



E 



us 



^ 4+Qa' 
ka/3 



, (2) 



where is the band dispersion, U the interaction 
strength, a 3 the third Pauli matrix, while the order pa- 
rameter S is defined by (I//V)( J2 ka 



P l 'k+Qa a a 



with 

N the number of sites. This one-body Hamiltonian can 



be diagonalized by the Bogoliubov transformation 

7ka = "kCka + Wk^(CT 3 ) Q /3Ck+Q/3 , 

P 

7ka = WkCka - Wk^(cr 3 )Q/3Ck+Q/3 ■ 

P 



(3) 



To avoid double counting, k is restricted to the magnetic 
zone. The superscripts c and v refer to the conduction 
and the valence bands split by the exchange Bragg scat- 
tering from the SDW. For simplicity, we assume perfect 
nesting e k = — £ k+Q- In this case, the coefficients of 
the transformation are given by u| = + ^k/^k)]. 



< =[!(!- ek/^k)], El = (e 2 + A 2 ), where A = US/2 



is the SDW gap parameter™. The diagonalized Hamilto- 
nian is given by H = £* Q #k(7 ka 7 ka " 7 ka 7L*) whcrc 
Y% means that the sum extends over the magnetic zone. 
The single-particle energy spectrum is given by ±£k and 
the SDW ground state for a half-filled band is defined by 
7icqI^) = 7kcJ^) = 0, which may be found by filling the 
vacuum with valence-band quasiparticles : 



l"> = Il^cL - £ 4 + Q^|J |0> 



(4) 



ka 



= ]J(i; k -u k £4 +Q/3 CT^ Q c kQ )4 a |0). (5) 

ka P 

The last form makes the analogies with the BCS ground 
state clear. For example, there exists Andreev-like reflec- 
tions at AF-N interfaces 14,15 . The ground state contains 
coherent particle-hole pairs. To clarify this, we perform 
a particle-hole transformation for states that are in the 
first magnetic Brillouin zone. Recall that destroying an 
electron in a state, creates a hole in the corresponding 
time-reversed state. For a spinor this can be achieved by 

Ckf -> ^tki and C ki - 

takes the form 



h_^. The ground state then 



\ty = I]> k - u kC k+Qt ^L k; )K - Wk4+cu /l -kt) 

k 

(6) 

where c k +Q Q \0) h = and /i kQ \0) h = 0. The particle- 
hole pair is in a triplet state with vanishing net spin pro- 
jection along the quantization axis, has no charge and has 
a wave vector equal to the antiferromagnetic wave vector. 
In the case of a FM, that wave vector would vanish. 

In closing this section, note that for a Neel vector ori- 
ented in some arbitrary direction with respect to the 
quantization axis, the variational wave function Eq.([5]) 
may be written 
* 

|fi> - Il( v k -«k£ 4+Q^4yf/t 7 c k7 ) |0) h • (7) 



P$7 



where 
U(8 



e -#/2 cog (0/ 2 ) _ e -i<t>/2 sin (0/ 2 ) \ 
e 4 */ 2 sin (0/2) e^ 2 cos(6>/2) J 



(8) 
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In a BCS superconductor, coherence of the pairs is re- 
flected by the fact that they all come with the same phase 
in the BCS state ]lk( M k-e # Wk4t ct -k4) 1°) - The u C 1 ) 
rotation e 1 *, identical for all pairs in this case, has its 
analog in the SU (2) rotation Uspcr^ = s ■ cr^with 
s = (sin#cos(/>, sin#siri(/>, cos0) that, in the antiferro- 
magnetic case, is applied to all particle-hole pairs in the 
same way. 

Although our results apply to any dimension, we often 
quote wave vectors for two dimensions to simplify the 
discussion. 



B. Tunneling Hamiltonian 

We choose to define the spin quantization axis along 
the direction of the instantaneous staggered magnetic 
moment S. Since we are interested in the case where 
the moments of the two AF are non-colinear (the analog 
of a phase difference), the spin quantization axis will dif- 
fer on each side of the junction. It is thus necessary to 
include a unitary transformation in spin space, defined 
by U(0, <fi) in Eq.([S]) above, to account for the fact that 
a spin up on one side of the junction is not the same as 
a spin up on the other side of the junction. The angles 
(9, 4>) correspond to the orientation of S# of the right AF 
expressed in the coordinate system of the left side of the 
junction; in Cartesian coordinates = |Si|(0,0, 1) and 
Sjj = | Sjj | (sin 9 cos 4>, sin 9 sin (j>, cos 9). The annihilation 
operators transform as follows 

^qcr U & a ' d,Q a i (9) 

where c? qCT / destroys a particle with spin quantization axis 
along the direction of the Neel vector on the right-hand 
side of the junction. The summation over repeated in- 
dices is implied. The tunneling Hamiltonian then reads 

H T = (1/A) (tkq4r I W*ia' + h.c). (10) 

The spin flip terms come purely from the choice of differ- 
ent quantization axes on the left and on the right. There 
is no real spin flip in the tunneling process. 

By choosing the quantization axis along the Neel vec- 
tor, we performed the analog of a gauge choice. The 
choice of a different basis (quantization axis) on either 
side is reflected completely in the XJ(6,<f>) appearing in 
the tunneling Hamiltonian. The angles in the unitary 
transformation U(0, <fi) are the non-Abelian analog of 
the phase in the ordinary Josephson effect. In the lat- 
ter case, electromagnetic fields coupled to charge affect 
the phase difference e - ie t/ hc through a gauge transforma- 
tion — > An + 9 M £ on one side of the junction. Since 
we can write U(#, </>) = e - jcr ^/ 2 e ~ lcr i/ e / 2 f or example, we 
see that in our case the gauge fields 9, (f> are coupled to 
spins, as is the magnetic field through the Zeeman effect. 
Such a field B yields a unitary time evolution e I£,AtBB crt / 2 



where where g is the gyromagnetic ratio, \x B is the Bohr 
magneton and t the time. Here it is the physical field and 
not the gauge field that affects the generators of SU (2) 
rotations. 



III. SPIN JOSEPHSON EFFECT 

To obtain the usual Josephson effect in superconduc- 
tivity, one computes the time evolution of the number 
operator, which is conjugate to the phase. The phase 
is the infinitesimal generator of the broken U (1) rota- 
tions. For a broken symmetry that is the staggered mag- 
netization along z, let us first consider the infinitesimal 
generator of SU (2) rotations along x, S x (q = 0) .Since 
the commutator [S x (q = 0) , S v (q = (n, n))] is propor- 
tional to the order parameter, a constant, this means 
that the analog of the number operator in this case is 
S v (Q — ( 7r 7 7r )) -By considering S v (q = 0) for infinites- 
imal rotations along y,one could have concluded that 
S x (q = (tt, 7t)) was the analog of the number operator—. 
Furthermore, S z will change directions. This means that 
one will find the analog of the Josephson effect, both DC 
and AC, by writing down the equations of motion for 
the vector S (q = (n, 7r)) .This is what we proceed to do 
in this section. 



A. Steady state 

The staggered magnetic moment operator in the left 
lead is S L = (h/2) J2k a fi ci +Qo ,<T a pc k p where <x is 
a vector of Pauli matrices. In the broken symme- 
try state, the time evolution of Sl due to Hl is 
negligible. Since we also have [§£ , Hr] — we 
find dS L /dt = (l/ih)[S L ,H T ], and thus dS L /dt = 

-(t/2-W) Ekq E Q( 35 *kq4 +Qa d q5 - h.C.) whose 

average Sl (t) = (dS l / dt) is given by 

kq ck/3<5 

(ii) 

where (...) is the thermal statistical average with the full 
density matrix. Performing first order perturbation the- 
ory using Ht as the perturbation, one obtains 

/t 
dt'([cl +Qa (t)d^(t),H T (t')}) , 

(12) 

where the average (...)o is computed with H — Hl + Hr 
(the unperturbed part of H ) . The operators on the right 
are in the interaction representation. 

Once the commutator is evaluated in Eq. (|12l) . 
one of the factorizations of the four-point correla- 
tion function involves products of correlation func- 
tions on the left and on the right leads such as 
(Ck + Q Q1 (*)cka(tO)o(<iq+Q«(*)^ 1 5(0)o- Such correlation 
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functions would vanish in a normal paramagnetic state. 
They are non-zero because of the broken symmetry. 
They represent interference in the tunneling process 
between momentum k + Q spin up particles and mo- 
mentum — k spin-down holes, in other words tunnel- 
ing of charge zero spin one S z — coherent particle- 
hole pairs that have finite momentum and are present 
in the ground-state Eq.©. In the case of the ordi- 
nary Josephson effect, one would find terms such as 

( c ia( t ) c -k-a( t '))o( d <iv( t ) d -<i-v( t '))o that represent tun- 
neling of coherent Cooper pairs. 

In order to compute the averages (...)o in the bro- 
ken symmetry states, we invert the Bogoliubov trans- 
formation Eq. (J3J) . Assuming i kq = tk q +Q = ^k+Q q = 
ik+Qq+Q, we find 

E^q4 +Q ^ = E E ^ q (rg q ) ii7 t + a 7^. (13) 

kq kq i,j£{c,v} 



where we defined 

( r kq)cc = (Uktiq + vlaVkUq + O^UkWq + ^la^SS^Vq) 
( r kq)c« = ("kUq - <7 3 aa U k Uq + Crf 5 U k U q - <4 a°««Wk W q) 
( r kq)fC = ("k-Uq + vlaVkVq ~ CT^U k U q - CT^CT^fcUq) 
(rg*) v „ EE (« k W q - CT^ Q U k W q - CT^t> k U q + ^CT^UkUq) . 

(14) 

Similarly, let Ht denote the part of Ht that does 
not commute with c k+ Q Q rf q 5. It can be written as 

Ht = (l/^Ekq^E^U^tkq^k^y^WL; Sub- 
stituting these expressions into Eq. (1121) and taking into 
account the fact that to this order in ikq the unitary 
transformation U (t') in the tunneling matrix element 
can be evaluated at t' = t since Hl and Hr do not change 
the quantization axis, one finds 



Si(t) 



1 * 

^EEE Im 



kq a/35 cr8' 



^|*kq| 2 ! Q. / 3U^U* ( 5/ 



dt , e -o+(t-t'). 



(15) 



S k > Q (*' 







where ^kfq) are ^he Keldysh Green func- 
tions in the left (right) lead. Their definitions 
are Sk(qw(M') *(7k + (q)^(* / )7t(q) a («)) and 

Sk(q),afl(*»*') = - z (TL(q)a( i )7 k t (q)^( i ')), respectively. 
Explicitly, 

££.'(*' ~ *) = - exp[-^|(t' - i)//i]<W , 

Sk^(*' ~ <) = exp[-^|(^ " *)A]<W . 

where / is the Fermi function. Eq. (fT5"j) for the staggered 
magnetic moment current through a tunnel junction is 
general. A bias could be included. We assume that there 
is no bias so there is no incoherent single-particle tun- 
neling across the antiferromagnetic gap. Integrating over 
t — t' and performing the spin sum in (1151) , one finds 

Sl = I c s L x s R , (16) 

where &l(r) = Sl(.r)/|S£(#)|.To obtain the correct sign, 
one must take into account in the Fourier transforms that 
the last site to the left and the first site to the right do 
not belong to the same sublattice. In the above equation, 
I c is defined by 

_ 8AiAfi A |2 /(£k)-/(-£k) 



with P indicating principal part. A similar expression is 
found for the equilibrium spin current in the case of fer- 
romagnetic tunnel junctions. Note that the sine function 
present in the standard Josephson case is replaced here 
by a cross product, which is a direct consequence of the 
vectorial nature of the order parameter. 

For a symmetrical junction (Ax, = Ar), the same as- 
sumptions and procedure as RefJ^ lead to the following 
analytical result 

Ic = ^R~ X A{T) tanh(f/3A(T)) , (17) 

where R — hj [ak^D 2 \t\ 2 ^j is the (zero-temperature) 

normal-state resistance of the junction with D the den- 
sity of state, which is assumed to be a constant. This 
expression for the temperature dependence of the criti- 
cal current has the same form as that obtained by Ambc- 
gaokar and Baratoff 1 ^ for a BCS superconductor, which 
is not surprising given the formal analogies^. 

By symmetry, the time derivative of the staggered 
magnetic moment on the right lead can be obtained by 
interchanging the L and R indices in Eq. (p~6|) . As a con- 
sequence, the staggered magnetic moments of the two 
AF precess about their (constant) sum S^, + Sr at a fre- 
quency w = I c \Sl + S_r|/|Sl||S_r|. 
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B. AC spin Josephson effect 

In the ordinary Josephson effect, the electromagnetic 
gauge potentials enter directly in the argument of the sine 
function. The present case is different. Each magnetic 
moment associated with a spin couples to the magnetic 
field through the Zeeman term (Hz = — ,g/i B B • S) where 
g is the gyromagnetic ratio and fi B the Bohr magneton 
(we neglect terms coming from orbital motion^). 

Considering magnetic fields Bl and H R applied re- 
spectively to the left- and right-hand sides of the junc- 
tion, the Heisenberg equations of motion lead to the fol- 
lowing equations of motion for the order parameters : 

Sl = -gfx B B L x Sl + I c §l x sr, 

Sr = -g^ B B R x S R + I c s R x s L . (18) 

The first term on the right side of the equality is merely 
the contribution of Hz to the Heisenberg equation of 
motion. The second term is the tunneling contribution 
and has exactly the same form as the one we have already 
computed in the zero- field case, namely Eq.([l6]). 

To prove the correctness of the 7 c sl x s r terms, we 
must return to the contribution from the first order in 
perturbation theory term Eg. (1121) . The expectation val- 
ues are now evaluated with the unperturbed Hamilto- 
nian Hl — <?/i s BL-S + H R — g{j, B H R S. Contrary to the 
case where no magnetic field was applied, the quan- 
tization axes now precess and one must worry about 
the time dependence of the rotation matrix U\, ,{t') in 

terms of the form $ g , (t')Ug, a , (t') c kcr' it') ■ That rotation 
matrix is no longer evaluated at time t. The creation- 
annihilation operators are in the interaction representa- 
tion and the original Hamiltonian commutes with the 
Zeeman Hamiltonian. Hence, it is possible to factorize 
the time evolution due to the magnetic fields to obtain 
QsS') = ^s(t')Afl(t') and c w (*') = A^(*')cL(*') 
where 

A^ L) it') = exp[i 9f i B B R{L) • at'/2] a , a (19) 

and the superscript on creation-annihilation operators 
indicates time evolution with Ho = Hl + H R . Since 
A^, a (f) = k L a , a it' - t) k L a - it) we can also write 

c k At') = ^,At'-t)A^(t)c^(t') (20) 

and define creation-annihilation operators whose spin ba- 
sis is rotated by the magnetic field, 

4 ff (t')^A^(f)4(t') 

so that 

c ka >(t') = A^(t'-t)cUt'). (21) 



The superscript t in (t') reminds us that the spin basis 
is that defined by the Neel vector at time t. That basis 
rotates because of the magnetic field. Somewhat redun- 
dantly, we also use a and a' to denote the spin basis 
on the left of the junction at time t and t' respectively. 
Physically, the last result tells us that the time evolution 
of the basis due to the magnetic field can be taken into 
account separately. 

We now need to work out the corresponding transfor- 
mation laws for the unitary transformation lit, ,(t'). This 
operator changes the basis from left (a 1 ) to right (5') at 

time t' . The only time dependence in lit, , (t r ) comes from 
that of the basis on the left and on the right due to the 
magnetic fields, it is not an additional time dependence. 
Hence, the transformation law, which we prove in more 
details below, is 

uLs') = A $s (*' - *) uimtl (*' - *) • ( 22 ) 

where we let 5 and 5' denote the spin basis on the right at 
times t and t'respectively. Given the transformation law 
for the operators on the left Eq. ([2"T|) and the analogous 
one for the operators on the right, we will thus find 

^At')Ul a ,it')c ka ,(t') = d%it')ulit)ci a (t') (23) 

i.e. all spin indices are at time t on the right-hand side of 
the equality. Therefore, it is possible to factor the unitary 
transformation Ul it) out of the integral over t' in the 
perturbation formula Eq. (TT2"| . The rest of the calculation 
then just follows the same path as in the zero-field case, 
leading to the above equations of motion Eq. (fT8|) . 

We now complete the proof on the transformation 
properties of the unitary matrices in Eq. (|22p . First, ex- 
press the left-hand side creation operators alternatively 
in the right and in the left-hand side basis. More specif- 
ically, left multiplying by W the change of basis Eq.©, 
the annihilation operators of left-hand side can be ex- 
pressed in the spin basis of the right-hand side 

ul a/ (t')c klT ,(t') = c k5 ,(t')- (24) 

Now we change to the right-hand side basis at time t with 
Eq.® 

Ul a , (tO cw(i') = A?, s (t' *) cUf ) (25) 

and we return to the creation operators in the spin basis 
appropriate for the left-hand side at time t. 

U\, a , it') cw (0 - A* it' - t) Ul (t) cL(t'). (26) 

We arrive at our final result by inserting the identity 
and changing the spin basis from time t to time t' , using 
Ea.(|2l]) again 
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Ut, a , (*') (f ) = Af, s (*' - *) Ul (t) (a£, (t> t) k L a , a (f t)) 
= (A& (f t) Ul (t) (f - *)) cw (i'). 



(27) 
(28) 



Comparing the left and right-hand sides, this completes 
the proof of the transformation property for the change 
of basis from left to right, Eq. ([23|) at different times. 

We close this section by solving and interpreting the so- 
lution of the equations of motion Eq. (fT8| . If B^ = B/v, = 
B, one finds that in the rotating frame u = {x',y',z'} 
defined by du/dt = — g/i B B x u, and Sr still pre- 
cess about their (constant) sum S = Sl + Sr at a fre- 
quency l}q. Returning to the static frame, Eq. (|18l) gives 
dS/dt = — ,g/i B B x S so that and Sr undergo a mo- 
tion of double precession. 

In the ordinary Josephson effect, a constant electric 
potential difference V leads to a phase difference <p that 
depends linearly on time (fp = 2eV/H). In the present 
case, an analog is found by computing the time depen- 
dence of the relative orientation i? between and Sr. 
Using Eq. (fT%|) to compute d(cos$)/cft = d(sL 4 SR)/dt, 
we obtain 



= -g/i B SB • e(t), 



(29) 



where <5B = Bl — B^ and e(t) = SR(t) x Si(t). As in 
the ordinary Josephson effect, this expression does not 
contain explicitly the tunneling matrix element. Con- 
trary to the ordinary Josephson effect, this equation is 
non-linear. Solving it numerically along with Eq. (|18[) 
in the case where the magnetic field vanishes on one side 
of the junction, one finds that the angle between and 
Sr behaves as a sine-like function of time. The presence 
of an additional constant magnetic field B throughout 
the system adds a beat to this sine-like behavior. The 
gyromagnetic ratio g is material dependent. No such 
non-universal constant appears in the ordinary Joseph- 
son effect. The above discussion of the AC effect can be 
transposed for FM by replacing the staggered magnetic 
moment by the uniform one. 



IV. EXPERIMENTAL DETECTION 



order of the conductance quantum, Rr 1 = 2e 2 /h, one 
finds a value of ojq of order 10 14 Hz, in other words in 
the visible. This frequency would be higher than the 
single particle gap and so would lead to much damping. 
Resistances that are many orders of magnitude larger are 
thus needed to bring its value down. 

To detect the spin Josephson effect even in the spin 
1/2 case, we suggest to add an additional tunnel junc- 
tion to a ferromagnet. In other words, consider an 
AF/I/AF/I/FM junction. Indeed, let us begin just with 
the AF/I/FM junction. The ferromagnetic material can 
be described by the Stoner Hamiltonian 



H 



C^-h o 
o c^ + h 



(30) 



where £ q = e q — \l and h is a molecular field. Using the 
same tunneling Hamiltonian as before and following the 
same procedure, we find that 



S L = I c Mr x S L 

Mr = I C S L X Mr 



(31) 
(32) 



where Mr is the average magnetization of the ferromag- 
net and 

f(^-h)-f (-Ik) 



C q - h + E k 



(h -*-h)\. (34) 



Since the uniform magnetization of the ferromagnet pro- 
cesses due to the Neel order parameter, it should be pos- 
sible to detect the precession of the FM moment by stan- 
dard magnetic resonance experiments. The addition of 
a second AF junction, as suggested above, should affect 
the ferromagnetic resonance. 



In principle, to detect the effects described above, 
one could use antiferromagnetic resonance 1 ^. This phe- 
nomenon is based on an anisotropy field resulting from 
the spin-orbit interaction. However the effect is absent 
for a spin 1/2. Assuming however that the same equa- 
tions would hold for S = 1, it becomes possible to detect 
the effect by antiferromagnetic resonance 1 ^. In the limit 
where Sl and S# are nearly colinear, the effect should 
therefore lead to a uniform shift of order ujq in the an- 
tiferromagnetic resonance frequencies of each AF. With 
= U = 2eV and a normal state conductance of the 



V. GENERALIZED JOSEPHSON EFFECTS 
BETWEEN SYSTEMS WHERE 
ANTIFERROMAGNETISM AND d-WAVE 
SUPERCONDUCTIVITY COEXIST 

In this section, we apply the same approach to study 
tunneling currents across a junction between two mate- 
rials that are d-wave superconducting antifcrromagnets, 
i.e. with homogeneous coexistence of antiferromagnetism 
and d-wave superconductivity. Since in the end we ex- 
pect symmetry considerations instead of details of the 
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Hamiltonian to be the dominant effect, we start with the 
simplest mean-field model that yields coexistence of the 
two order parameters. Then we show that for coexisting 
orders, the ordinary Josephson effect can be modulated 
by the antifcrromagnetism and that the spin Josephson 
effect can be modulated by the superconductivity. These 
effects occur already to leading order in the square of the 
tunneling matrix element. 



different interactions may be defined in the following way: 



c„: c it 



( n i<r) = 2 + °" mc °s(Q ■ n) , (39) 



(40) 



9 ^29(8)(ci+6,tCi,i + Ci+8,iCirf) = t cos(Q • Y t ) , (41) 



A. Model 

To take the simplest possible case that has the appro- 
priate symmetries, we model the AF/SC coexistence on 
each side of the junction by the following phenomenolog- 
ical mean-field Hamiltonian^ : 

- vj2(&y^) + k - c -) - ^£( a u a m) + h.c.) 

i i 

(35) 

The possibility of a spin-triplet pair amplitude (A^j) 
must also be considered in general for self-consistency2i. 
The singlet A^.i and triplet A t .i pair destruction opera- 
tors are defined as 

A <« = 2^ 9 ^ S ^ Ci+s ^ Ci ^ ~ c »+ 5 .4- c «.t)i (36) 
s 

A M = \ 9( s )( c i+s,fCi^ + c i+s ^c l:t ) (37) 



and the structure factor g(S) is chosen to have a d- wave- 
like form such as 



1/2 if 5 = (±1,0) , 
-1/2 if 5 = (Q,±1), 
if otherwise . 



(38) 



The strength of AF, SC, and spin-triplet pair interac- 
tions is governed by U, V, and W, respectively, which 
are all positive in our study. £ k is a tight binding en- 
ergy dispersion given as £ k = —2t(cosk x + cosk y ) — 
At' cos k x cos k y , where t and t' are hopping constants for 
nearest neighbors and next nearest neighbors, and fi is 
the chemical potential, which controls the particle den- 
sity n. 

For completeness, we recall the results of Ref. [2(3 that 
we need. In the mean-field approximation, three different 
order parameters m, s, and t corresponding to the three 



where Q is the (commensurate) antifcrromagnctic wave 
vector equal to (tt, tt) in two dimensions. Now the mean- 
field Hamiltonian Hmf is quadratic in the original elec- 
tron operator and in terms of a new four component field 
operator ?/> k it becomes bilinear 



where 



H C M = Wc M kV>k + #0 , 



V'k = ( c kt> c -ki> c k+Qt> c -k-Q±) 



(42) 



(43) 



The matrix M k is given as 



where 



/ e k Vs<f>(k) 

Vs<j>(k) -e k 

-Urn -Wt<j>{k) 

\Wt(/)(k) -Urn 



-Urn 
-Wt<f)(k) 

e k+Q 

-Vs(j>{k) 



ek = Ck - V 



Wtcj>{k) \ 
-Urn I 

-Vs(j){k) I ' 
-«k+Q / 

(44) 



(45) 



and (f>(k) = cosk x — cosk y is the Fourier transform of 
g(S). The constant energy shift Eq depending on to, s, t 
is given as 



E = N{Um 2 + Vs 2 + Wt 2 - /j) , 



(46) 



where N is the total number of lattice sites. The energy 
eigenvalues of M k yield four energy dispersions ±E±(k) 
(with E±(k) > 0) l2(fcl 

E ± (k) ={(4 + 6 k+Q )/2 + (Urn) 2 + [Vs^k)} 2 

+ ^(k)fi S (k)} 1 / 2 , (47) 

where g(k) is given as 

g(k) ={(4 - £ k+Q ) 2 /4+ (6k - £k+Q ) 2 [^(k)] 2 

+ ((ek + £k +Q )(t/m) + 2{Vs<f ) (k)}{Wtcf>(k)}) 2 } 1 / 2 . 

(48) 

When s = t = or m = t = the energy eigenvalues 
reduce to those obtained from the SDW or BCS factor- 
izations of the corresponding interactions. 
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The transformation matrix A relating the original c 
operator to the eigenoperator 7 is defined by 

^k = X>y(k)7 ik . (49) 

3 

All the elements of the matrix A are given in the appendix 
of Ref. [20|- The columns of the matrix A are the eigen- 
vectors of the mean-held Hamiltonian with the corre- 
sponding eigenenergics ordered as follows, = +E + (k), 
El = -E + (k), El = -J5_(k), Ei = +£L(k). The vac- 
uum state of "0k i s obtained by filling the vacuum of the 
original operators Cko- with down spin electrons. As usual 
in the grand canonical ensemble, the ground state is then 
obtained by filling the vacuum of ip k with all the nega- 

I 



tive energy quasiparticles i.e. with j 2]i and 7g k whose 
energies, = —E + (k) and El = — £L(k) respectively, 
are negative for all k. One recovers the BCS and SDW 
ground states in the appropriate limits. In general that 
ground state contains triplet pairs. 



B. Combined spin Josephson effect and ordinary 
Josephson effect in superconducting 
antiferromagnets 

We now turn to the calculation of the spin and charge 
tunneling current across the junction. In terms of the 
7's, the equation for dS^/dt can be rewritten as 



sow = jvEEEM^p^/ 2 ^^^^ 

k,q i,j (50) 

+ ^U^f^ t (7ik(*)7iq(*)> +^U W f« u ( 7ik (t) 7 J q (t))}] , 



where the phase difference and the difference in quanti- 
zation axis have been gauged into the tunneling matrix 
element and where we defined 

fkVr = t^( k ) + 4K(k)HA/(q) + ^3i(q)] , 

= [^( k ) + 4w(k)][4y(q) + ^-(q)] . 

f LVt = ^i(k) + A 4i (k)]L4 Xi (q) + A 33 (a)] , 
= [4«(k) + A 4l (k)][A 2j (q) + Ay(q)] . 

We compute each average in (|50|) to first order in the 
tunneling Hamiltonian as in Eq. (|12[) which r equires Ht 
to be also rewritten in terms of the 7's. However, un- 
like the previous case of itinerant antiferromagnetism, 



all the terms in Ht must now be included since uncon- 
ventional pairing is allowed by the presence of supercon- 
ductivity. There will appear two-point correlation func- 
tions of the form (cj co .(t)cl k _ (T (i / )d q(7 (t)d_ q _ .(*')}o) cor- 
responding to the coherent tunneling of Cooper pairs, 
and these will acquire a factor exp^Ayj), where Aip is 
the phase difference between the superconducting order 
parameters of the two leads. Taking this into account, 
the calculation of the staggered magnetic moment cur- 
rent leads to 

8 G (t) = (l c + Jc cos A(p)s D x § G , (52) 

where 



1 * 



k,q 



F? 
^q 



El 



1-/(4) -/(££)' 



El 



(53) 



and 



E 772 Hl^ql 2 
k,q 



-pV TV T3 
1 kqtt 1 kqU r 



f(K) - f( E l) 



F? 
^q 



El 



i-fm)-f(Ei) 



El 



(54) 



with defined at the end of the previous subsection. 

Similarly, we obtain the charge supercurrcnt (or 
Josephson current) across the junction by computing the 
statistical average of the time derivative of the number 
operator on the left lead, Nq — Sko- c iko- Ckcr: 



N G (t) = (l c + J c cos 8) sinAip 



(55) 
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with 



'N 2 h ^ 

k,q 



rkq| J- kqft kqU 



/(^q) ~ f(K) 



F? 

^q 



1 - /(££) - /(£«) 



£q+^ 



(56) 



and J c is obtained from I c by changing the sign be- 
tween the two quotients in the brackets. Note that 
J c = 2J c /h.We discuss the consequences of this equal- 
ity in the next section. 



VI. DISCUSSION AND CONCLUSION 

The robustness of the spin-Josephson effect is related 
to the rigidity of a broken symmetry, in analogy with 
the ordinary Josephson effect. From a phenomenological 
point of view, the common thread between the gener- 
alized Josephson effects described in this work is that 
the effective Ginzburg-Landau free energy, or the ef- 
fective action, contains interaction terms proportional 
to product of the order parameters on the left and on 
the right of the tunnel junction. Symmetry imposes 
(^l^Jj + h.c.) oc cos(A</>) in the superconducting case 
and Sl (Q)-Sr.{— Q)in the antiferromagnetic case, with 
Sl (Q) the Neel vector. The latter case is completely 
analogous to the ferromagnetic one where one finds the 
scalar product of the uniform magnetizations 4-6 . The 
corresponding equations of motion involve sin (A^) in 
the superconducting case and its natural generalization 
in the vector case, namely the cross product. The result- 
ing antiferromagnetic spin current leads to a precession 
of the order parameters about their sum on either side 
of the junction. In the case of a spin 1/2, this precession 
can be detected experimentally by comparing the mag- 
netic resonance of the ferromagnet in an AF/I/FM junc- 
tion with the magnetic resonance of a AF/I/ AF/I/FM 
junction, as we have shown. We speculate that similar 
effects exist for spin one antiferromagnetic junctions in 
which case the effect could be directly observed through 
antiferromagnetic resonance^. 

In the mean-field solution, the direction of the triplet 
in spin space is locked to that of the antiferromag- 
netic order parameter. For a d -wave superconduct- 
ing antiferromagnet, there is an additional cross term 
in the Ginzburg-Landau free energy that comes from 
the triplet component W*of the order parameter, namely 
( W £ • + h.c.^j that leads to a) a modulation of the 

critical charge current by the dot product of the antifer- 
romagnetic order parameters and b) a modulation of the 
critical spin current by the superconducting phase differ- 
ence. An analogous effect was found for superconducting 
ferromagnets^. The equality J c = 2J c /h found at the 
end of the previous section is a direct consequence of the 
fact that both cross terms come from a single term in 



the effective Ginzburg-Landau free energy. There is thus 
a reciprocity between the modulation of the Josephson 
charge current by the antiferromagnetic order parameter 
and the modulation of the antiferromagnetic Josephson 
spin current by the superconducting order parameter. 

An additional analogy between the ordinary Joseph- 
son effect and the spin-Josephson effect resides in gauge 
transformations. In the ordinary Josephson effect, a 
gauge transformation allows the phase difference to be 
reported entirely on the tunneling matrix element. Simi- 
larly, in the spin Josephson effect, one can make a choice 
of quantization axis (analogous to a gauge choice) that 
transfers to the tunneling matrix element everything that 
is related to the difference in order parameter orientation 
across the junction. 

We learned from the microscopic calculation that tun- 
neling of Cooper pairs is replaced by tunneling of the cor- 
responding condensed objects in the antiferromagnetic 
case, namely spin-one, charge zero, particle-hole pairs, 
with zero projection along the direction of the antifer- 
romagnetic order parameter. Alternatively, and also in 
analogy with the superconducting case, the terms that 
lead to the spin-Josephson effect represent interference 
in the tunneling process between momentum k + Q spin 
up particles and momentum — k spin-down hole. In su- 
perconducting antiferromagnets, both types of condensed 
objects, Cooper pairs and finite Q particle-hole pairs, 
are present hence one sees both Josephson and spin- 
Josephson effects. In addition, in superconducting anti- 
ferromagnets a pure singlet state between nearest neigh- 
bors is impossible because of the broken spin symme- 
try between the two sublattices. Hence there are always 
triplet Cooper pairs with total wave vector Q21. It is 
the tunneling of this type of Cooper pair that leads to 
the cross term between charge and spin-Josephson effects 
mentioned above. Indeed, a triplet Cooper pair carries 
both phase and spin information and hence its tunneling 
can influence both spin and charge and can be manipu- 
lated by fields that couple to either spin or charge. 

One of the most far reaching consequences of the 
Josephson effect in ordinary superconductors is the AC 
effect where phase difference between the superconduc- 
tors increases linearly with time due to an applied voltage 
difference. The corresponding frequency of the current is 
proportional to the voltage drop through a material inde- 
pendent universal constant 2e/h. In the case of the spin 
Josephson effect, even without external field there is a 
precession frequency. In the presence of a uniform mag- 
netic field, the two Neel vectors precess around their sum 
while that sum precesses around the magnetic field. The 
correct analog of the AC effect corresponds to adding 
fields that will lead to a time variation of the angle be- 
tween the the two order parameters, the analog of the 
phase difference. Adding a uniform magnetic field does 
not achieve this. It is only in the presence of a mag- 
netic field difference between the two sides of the junction 
that the angle becomes time dependent, as described by 
Eq.(gH]), i9(t) = -g^ B 5B-e(t) with 5B the magnetic field 
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difference and e(t) = Su(t) x &L(t)- If the correct order 
parameters are used, this result is valid for tunneling be- 
tween ferromagnets as well. As far as we know, no analog 
if the equation = — g/i s (5B • e(t) has appeared in the 
literature on spin-Josephson effect. Contrary to the or- 
dinary Josephson effect, the proportionality constant in 
front of is not universal since g is material depen- 
dent. This is because the coupling of the magnetic field 
to matter is not a pure gauge coupling. In the presence 
of spin-orbit interactions, both charge and spin currents 
couple to the field. 

One possible application of our work is to use it to 
differentiate between homogeneous coexistence of anti- 
ferromagnetism and d-wave superconductivity and phase 
separation between these two types of order. It would 
probably be easiest to measure time dependent modula- 
tion of the " charge" Josephson critical current induced by 
a junction where a magnetic field gradient small enough 
not to destroy superconductivity is applied to the tun- 
nel junction. In the case of phase separation, the cross 
terms that we described are not present since they come 
from tunneling of a finite momentum triplet Cooper pair 
that exists only when there is homogeneous coexistence. 
The issue of homogeneous coexistence is important for 
the cuprates, the heavy fermions, the organics and the 



pnictides where antiferromagnetism and d-wave super- 
conductivity often seem to overlap in the phase diagram. 

Further possible generalizations of our work include 
a) Studying tunneling between states with more compli- 
cated order parameters b) Calculating the analogs of the 
effects that we found in the case where the tunnel junc- 
tion is replaced by a normal metal. It is known^r— that 
Andreev states in the superconducting case profoundly 
change the dependence of the current on the phase differ- 
ence, but not its periodicity. Since the analog of Andreev 
states exist in the antiferromagneti 4 , similar effects are 
expected for the spin-Josephson effect. 
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